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Abstract 

The aging dynamics of a simple model glass is numerically investi- 
gated observing how it takes place in the potential energy landscape 
V. Partitioning the landscape in basins of minima of |Vyp, we are 
able to elucidate some interesting topological properties of the aging pro- 
cess. The main result is the characterization of the long time behav- 
ior as a jump dynamics between basins of attraction of minima. More- 
over we extract some information about the landscape itself, determin- 
ing quantitatively few parameters describing it, such as the mean en- 
ergy barrier value and the mean square distance between adjacent minima. 

The specific way a glass or a supercooled liquid explores its potential energy landscape at 
different temperatures and in different dynamical regimes is of great importance to under- 
stand many properties of these systems (Stillinger, 1995; Debenedetti and Stillinger, 2001). 
Only recently the investigation of the multidimensional Potential Energy Surface (PES) as a 
function of particle coordinates has become computational achievable, and many numerical 
works have been devoted to study the PES of model liquids and glasses. The relationship 
between the topological properties of the landscape and the way the system moves in it, con- 
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stitutes a very promising viewpoint to better understand the thermodynamics and dynamics 
of supercooled hquids and glassy systems (Kurchan and Laloux, 1996). 

Only few studies have investigated the aging dynamics from the energy landscape point 
of view (Angelani et al. 2001; Kob et al. 2000; Mossa et al. 2001). The study of the 
off-equilibrium dynamics in glassy systems, and in particular in spin glasses, has allowed 
to obtain important information on the phase-space properties of the systems itself and 
to formulate off-equilibrium relations generalizing the usual equilibrium ones - generalized 
fluctuation-dissipation relation (Cughandolo, and Kurchan 1993; Bouchaud et al. 1997) - 
with the introduction of an effective temperature (Cughandolo et al. 1997). The investi- 
gation of the off-equilibrium properties seems to indicate the correctness of the conjecture 
stating the similarity of structural glasses with some spin glass model - spin glasses with 
one step replica symmetry breaking (Kirkpatrick and Thirumalai, 1987; Kirkpatrick and 
Wolynes, 1987; Pranz and Parisi, 1997). Recently the numerical investigation of aging in 
structural glasses has allowed us to bridge these two different disordered systems (Parisi 
1997; Barrat and Kob 1999; Di Leonardo et al. 2000; Sciortino and Tartaglia 2001;). Here 
we report a detailed analysis of the aging in the energy landscape for a simple model glass, 
a monatomic Lennard- Jones system modified in such a way as to avoid crystallization. 

The way the off-equilibrium dynamics takes place in the potential energy landscape V 
(viewed as a function of the 3N coordinates of the N particles) is investigated by mapping 
the true dynamics in the configuration space on to the "slow" points of V, i.e. the points 
in which |Vyp is minimum and the forces are reduced by several orders of magnitude (we 
call them saddles hereafter, even if, besids true saddles, they include also inflection points 
- Doye and Wales 2001; Angelani et al, in preparation). The basins of attraction of these 
points constitute a partition of the configuration space. This partition of the landscape is 
an extension of the partition in basins of minima proposed by Stillinger and Weber (1982, 
1984) many years ago, including in our case not only the minima but also the saddles of 
any order. The new partition has been very useful to elucidate many interesting features 
of supercooled simple liquids, the main result being the identification of the mode-coupling 
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temperature Tmct (Gotze, 1999) as the temperature at which the system has a crossover 
between a dynamics between basins of saddles for T > Tmct to a dynamics between basins 
of minima for T < Tmct (Angelani et al. 2000; Broderix et al. 2000; Scala et al. 2001). 
We use this new partition in the analysis of the off-equihbrium dynamics of a monatomic 
Lennard- Jones system with N — 256 particles. The particles interact via the usual 6 — 12 
Lennard- Jones potential modified in such a way as to avoid the crystalization occuring in 
standard LJ: V — Vlj-\-5V , where 5V is a small many-body term that inhibits cristalization 
(it depends on the static structure factor, see Di Leonardo et al. (2000) for details). The 
relevant quantitites are determined during the time evolution of the system after a sudden 
crunch (density jump at fixed temperature) from liquid to glassy phase (from po — 0.95 to 
Pi — 1.24 at fixed Tq = 0.5, where the glass transition temperature is Tg ~ 1.4 at the final 
density pi - all the quantities are expressed in standard Lennard- Jones unit). The evolution 
of the system is generated up to time t = 3 • 10^ by means of standard isothermal MD with 
time step r = 0.01. The measured quantities have been averaged over 80 samples. For each 
MD trajectory we have quenched a set of logarithmically equally spaced configurations to 
the corresponding minima and saddles following the steepest descent path respectively on 
the potential energy surface V in the former case and on the auxiliary potential W — \\W\^ 
(Weber and Stillinger 1985) in the latter. In this way, during the time evolution, we are able 
to follow the behavior of the sampled relevant points of the landscape: saddles and minima 
(besides the usual instantaneous configurations). 

In order to elucidate how the system explores the landscape during the off-equilibrium 
dynamics, we have first focused on the energy values. We have measured the potential ener- 
gies as a function of time (times are measured from crunch) along the dynamical trajectory: 
the energy of the instantaneous configurations, the energy of the saddles and the energy 
of the minima. The data are reported in Fig. 1 in a logarithmic scale. It is evident that 
the instantaneous and saddle energies have a rather sharp transition at time tc ~ 20. The 
meaning of this crossover will be clear soon, analyzing more carefully another topological 
quantity related to the landscape, the order of sampled saddles. We only notice here that 
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for times less than tc the instantaneous energy follows the energy of saddles, and the energy 
of minima remains nearly constant in this region (we note this fiat time-region corresponds 
to the constant energy region observed at high temperature in the analysis of minima in 
equilibrated liquids), while for times below tc all the energies (instantaneous, saddles and 
minima) have a similar almost logarithmic behavior, indicating that here the time evolution 
is driven by the inherent dynamics between minima: e(^) oc 7logio(t), with 7 ~ —0.025. 

Another remarkable point that can be observed in Fig. 1 is the similarity between the 
energy of the saddles and the difference between the instantaneous energy and the constant 
term ^Tq. This seems to indicate that the non-diffusive degrees of freedom give an harmonic 
contribution over the energy of saddles, slowly approaching their equilibrium value. 

A better understanding of the processes that are taking place is obtained analyzing the 
number of negative directions in the landscape near the sampled saddles. We define the 
saddle order Us as the number of negative curvatures at the saddle point, i.e. the fraction 
of negative eigenvalues of the Hessian of V. The time evolution of Ug during aging (Fig. 2) 
allows us to clarify the meaning of the crossover time tc'- it marks the time at which the 
system begins to visit extensively the minima, and the mean order ng falls below the value 
rig = 1. 

In the analysis of equilibrium properties of supercooled liquids, the quantity Ug vs. tem- 
perature allowed us to give a topological interpretation of the mode-coupling temperature 
as the temperature at which a transition takes place, from a dynamics between saddles to 
a dynamics between minima (Angelani et al. 2000). In a similar way, in the present study 
of the out-of-equilibrium properties of a simple glass, the same quantity Ug, now plotted vs. 
time t, allows us to interpret the time tc as marking the transition between the same two 
regimes: at short times, below tc, the system moves inside basins of attraction of saddles 
with mean order ng{t) > 1 decreasing in time, while at longer times, t > tc, the system starts 
to visit basins of minima traveling among them through low order saddles. The fact that 
the mean order Ug doesn't vanish for times t > tc, explains the observed energy difference 
between the energy of saddles and that of minima. 
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Analyzing how the system explore its potential energy landscape during the aging dy- 
namics we can try to infer some topological characteristic of the landscape. Matching the 
information about the order and energy of saddles and the energy of minima, we are able to 
evidence the possible existence of a typical energy barrier in the landscape. As reported in 
a previous work investigating the equilibrium supercooled dynamics (Angelani et al. 2000), 
also in the aging regime we discover a very simple topology of the landscape: saddles of 
order Ug lie above the underlying minimum by an energy elevation proportional to n^: 

es{ns) - = • . (1) 

In Fig. 3 the left hand side of Eq. 1 is reported as a function of saddle order Ug. The 
open circles are the data obtained from previous figures, in which the minima are defined 
by a steepest descent path from the instantaneous configuration. However, using a dif- 
ferent definition of underlying minima below a saddle (by quenching the saddle instead of 
instantaneous configuration - see the text below) one observes a very similar behaviour (full 
diamonds in Fig. 3) . This finding supports the picture of a very simple landscape with only 
a single energy barrier parameter AE which represents the energy gap between a saddle of 
order Ug and a saddle of order -|- 1: 

~ 10 . (2) 

Similar results were found in the supercooled equilibrium case, where a value of AE ~ 3.6 
was obtained for a less dense system p = 1 (in our case p = 1.24). 

Another interesting topological information about the landscape is the valuation of the 
mean distance between adjacent minima. In order to obtain this information we have 
determined, for each saddle point sampled during the dynamics, two underlying minima 
(here the term underlying means minima reached by a steepest descent path from a saddle 
configuration, weakly perturbed in order to escape from the instable point), {rj^^} and 
{r^^"*}, reached by quenching the saddle. In this case the above definition of minima, as those 
reached by saddles, seems to be the more appropiate, due to the fact that the minima sampled 
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during the dynamics starting from instantaneous configurations could be not directly related 
to the sampled saddles. Instead a direct quench from saddles assures the close relation in 
configuration space between the two points (we note that in the valuation of energies the two 
method lead to the same results, indicating that the energies are less sensitive to landscape 
details). The quantity d'^ defined as 



is reported in Fig. 4 as a function of the saddle order Ug. It is evident a linear relation 
(hne in Fig. 4) between the two quantities, suggesting that the descent path from a saddle 
of order towards the underlying minima can be interpreted as a random walk among 
saddles: at each step the order of the saddle decreases by 1 and the mean square distance 
increases by a constant value d^ (this characteristic, together with the single barrier energy 
parameter AE, is incorporated in the simple trigonometric model as defined by Madam and 
Keyes (1993)). The value obtained for the typical square distance between minima is: 



We note that this finding and the previous one concerning the mean barrier value (see 
Eq. 2), although obtained in the specific aging regime, are supposed to be properties of 
the landscape itself, independent on the different way the system explores it. Indeed the 
same characteristics were obtained by Angelani et al. (2000) who analysed the equilibrium 
dynamics in the supercooled regime (a quantitative match between the obtained values is 
not feasible, as the densities of the systems are different in the two cases). 

In conclusion, the topological investigation of the energy landscape of a simple model 
glass during the aging dynamics has allowed us to clarify many interesting properties, shortly 
resumed here after: 

• after a fast short time behavior, characterized by a dynamics between saddles, the 
system starts to explore minima decreasing logarithmically its energy (this latter is 
the interesting slow aging dynamics); 




(3) 




0.02 . 



(4) 
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• the instantaneous energies lie above the saddle energies by the harmonic term |T, 
indicating that the fast harmonic degrees of freedom equilibrate over the slow diffusive 
coordinates (related to the negative directions of the sampled saddles); 

• the topology of the landscape seems to be highly regular, with only a single energy 
barrier parameter AE ~ 10 and a single square distance between adjacent minima 
dl ~ 0.02. 
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FIGURES 




FIG. 1. Energies (per particle) during the aging dynamics. Instantaneous potential energy 
ei subtracted by the constant term 3/2T (open circles), energy of minima Cm (open squares) and 
energy of saddles (full circles). The solid line fits the long time data, with slope 7 ~ —0.025. 




FIG. 2. Time behavior of the mean order n, of the saddles. 
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FIG. 3. The energy elevation of saddles from the underlying minima as a function of saddle 
order n^. The open circles are the data obtained from Fig.l and Fig. 2 (the minima are defined 
starting from instantaneous configurations), while the open diamonds are obtained determining the 
underlying minima starting from saddle configurations. The fit line determines the energy barrier 
(expressed as total energy) AE = 10. 




FIG. 4. Mean square distance between adjacent minima vs. saddles order n^. The line is a 
linear fit, with slope = 0.02. 
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